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THE MODUIATION OF LOW ENERGY GALACTIC COSMIC 
RAYS OVER SOLAR MAXII^JM (CYCLE 20) 

M. A. I. Van Hollebeke, J. R. Wang, and F. B. McDonald 
NASA/Goddard Space Flight Center, Greenbelt, Maryland 20771 

This note is a preliminary report on detailed observa- 
tions of low energy (~ 6o MeV/nucleon) galactic cosmic rays 
made during the recent period of solar maximum. For this 
particular cycle a significant time lag between Intensity 
changes in the low energy and high energy components has been 
observed. There are a number of possible explanations of 
this "hysteresis effect". In the simplest form these include 
either changes in the form of the magnetic field power 
spectr\un or changes in the size of the effective modulating 
region, or both. 

The time histories of the 60 MeV/nuc galactic protons 
and alpha particles and of the high energy component as 
represented by the Deep River Neutron Monitor rates, for the 
period June I967 - October 1971, are shown in Fig. 1. The low 
energy data was obtained from the GSFC cosmic ray experiments 
flown on the IMP IV and V satellites. Special care was taken 
to Insure there was no solar particle contamination in the 
40-80 MeV/nuc energy interval. This was done by examining 
the spectra over the complete 4-80 MeV interval and selecting 
only those periods where the 4-30 MeV flux was low and 
characterized by a steeply falling energy spectrum. 

Examination of Fig. 1 shows that the low energy compon- 
ents and the neutron monitor rates display the same trends 
from 1967 to mid-1970. However, by early 1971 the Deep River 
Neutron Monitor rate has returned to its I967 level while 
the 60 MeV/nuc proton and alpha particle Intensities have just 
begun to show significant Increases. A similar behavior has 
previously been observed for 0.02 - 20 GeV electrons (Meyer 
et al 1971) and for 1-10 MeV secondaries neutrons (Verschell 
et al 1971). / 
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This "hysteresis effect" is more clearly defined by 
plotting the 6o MeV/nucleon protons and alphas versus the 
Deep River Neutron Monitor counting rate on a semi-logarithmic 
scale (Fig. 2a and h). It is seen that the proton data can 
he fitted hy two essentially straight and parallel lines. The 
top line covers the period from mid -1965 to late I 969 and 
the lower parallel line covers the period July 1970 to Sept- 
ember 1971 . The alphas (Fig. 2b) display an almost identical 
behavior except that the displacement of the two parallel lines 
at a given Deep River rate corresponds to a factor of 2 change 
in alpha flux compared with a factor of 2.6 for the protons 
(Fig. 2a). The period from late 1969 to July 1970 is one of 
transition with the data suggesting that the alphas "lead" 
the protons. Earl and Rygg (private communication) observed 
only a factor of "" 1.04 for protons above 750 MeV. These 
results show that there is a well defined hysteresis in solar 
cycle 20 and that this effect is strongly rigidity dependent. 
They also confirm the preliminary reports given independently 
by Van Hollebeke et al. (1972), Garcla-Munoz et al. (1972), 
and Lockwood et al. (1972). 

The question then arises as to whether the phenomena we 
report here is a characteristic feature of each solar cycle 
or is unique to cycle 20. For the maximum of the solar cycle 
19, Simpson (I963) reported a significant hysteresis during 
195'^-1962 between neutron monitor rates at Chicago and Climax 
(vertical cutoff rigidity of ~ 1.3 GV and ~ 3 GV respectively) 
and at Huancayo (vertical cutoff rigidity ~ 13 GV ) . However , 
this was reduced to almost negligible proportions when correc- 
tions due to the drift in the Huancayo aneroid pressure gauge 
were made (Simpson and Wang, 1970) . From a regression plot 
of the low-energy data obtained from polar^ high altitude 
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obs ervatlons (Neher, 1971 )j and Mt. Washington neutron monitor 
rates 5 Webber (I967) found no evidence for hysteresis over 
cycle 19, Thus in our sample of 2 solar maxima, we find one 
displaying a well defined hysteresis extending up to energies 
of several hundred MeV and an apparent absence in the previous 
cycle. In this regard, it is interesting to note that during 
cycle 19, the sunspot niunber peaked at a value of ~ 200 for 
only 6 months, while during cycle 20, it stayed at a maximum 
value of 105 for almost 3 years. 

On the other hand, by comparing the data from Geiger 
counter arrays on IMPs I, II, and III (protons > 60 MeV, 
electrons > 3 MeV) and the Deep River Neutron Monitor rates 
for the 1963-1966 period, Bala Subrahmanyan et al. (I968) found 
a well defined hysteresis effect over the last solar mlnlm\Am . 
These results were confirmed by Kane and Winckler (I969) using 
ion chamber measurements (protons and 'alphas > 12 MeV/nuc, 
electrons > 0.6 MeV) on OGOs I and III satellites. However, 
in a compilation of differential spectrum measurements from 
several different experiments, OGallagher (1969) discussed 
the previous results and report a negligeable effect over the 
same period. A further analysis of this part of the solar 
cycle should be made to resolve this discrepancy. 

The shapes of the low-energy cosmic proton and alpha 
spectra measured at solar minimvim have been accounted for by 
Goldstein et al. (1970) using a numerical solution (Fisk, 1971) 
to the spherically-symmetrlc Fokker-Planck equation which 
allows for the effects of convection, diffusion and energy 
loss. Their results are consistent with an interstellar 
differential spectrum proportional to (To + T)~^'^^, where 
To is the particle rest mass energy and T the kinetic energy. 
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The Fokker-Planck equation is given by { Parker ^ 19 ^^; Gleeson 
and Axfordj 1967) 
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^(e^P^r^t) = particle diffusion coefficient^ 

3 = particle speed in units of c 

P = particle magnetic rigidity 
We find that to theoretically reproduce from the above 
equation the observed proton and. alpha differential energy 
spectra, as well as the hysteresis effect over a solar cycle, 
it appears necessary to change either the rigidity dependence 
of k: ^ or the size of the modulation region as a function of 
time or both of these. In Table I we give two possible forms 
of K . These differ slightly from the one used by Goldstein 
et al. (1970). The form Kq- constant in r in the inner solar 
system has recently been advocated by Fisk et al. (1972). 

The values shown in Table I for A{t) and T(t) for case 
I, and m(t) for case II, are used with the interstellar proton 
and alpha differential energy spectra of Goldstein et al. 
(1970) to produce the dashed curves in Figs. 2a and 2b. These 
. parameters are changed smoothly over the indicated range for 
a given period. We emphasize that neither of these cases or 
the choice of parameters are unique and that both cases can be 
made to produce identical results by the proper choice of 
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parameterSo For example^ in case I we could vary D instead 
of A, Furthermore^ we could also consider a comhination of 
cases I and II. Other possible ways of producing a hysteresis 
effect by changing the rigidity dependence of K have been 
pointed out by Nagashima et al. (I 966 ) and by Urch and 
Gleeson ( 1972 ). 

It is interesting to note that in case I, y(t) changes 
only from 1.0 to 1.2. According to Jokipii (I 966 ), this 
corresponds to a decrease on the order of 20 ^ in the index 
of the magnetic field power spectrum. This would appear to 
be more physically reasonable than the large changes required 
in D(t) in case II. However the suggestion that the alpha 
particles lead the proton in the transition period^ making 
the lag time shorter for high rigidity, seems to favor a 
model depending on the size of the modulating region. The 
radial cosmic-ray intensity gradients predicted by these two 
cases are somewhat different. For example, the intensity 
gradient of " 60 MeV protons during late 1971 for case I 
and II are ~ 95^/A.U. and "" 62 ^/A.U. respectively. 

The fact that such a well defined hysteresis is observed 
over solar maximum in cycle 20 implies changes in the state 
of the Interplanetary medium. Magnetic field power spectra 
over this interval are not yet available. Examination of 
some parameters such as the solar wind speed and plasma ion 
density have been made over most of this period.it has been 
shown (Gosling et al. 1971) that the average plasma speed V 
remains approximately constant in time around a value of 400 
km/sec. However, the 27 day averages plot. of V shows a slight 
maximiun of "" 20 ^ above this mean value during the period from 
late 1967-1968 (Montgomery and Bame, 1972). As it can be 
seen in Figure 2 this small change during this period cannot 
produce the observed hysteresis but it could be one of the 
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factors affecting the power spectr-um of magnetic field 
fluctuation making (Table I) smaller. Further, since 
the condition for hydromagnetic instabilities (which presum- 
ably generate magnetic field irregularities) to occur depends 
strongly on the ion plasma density (Burlaga, 1971 , Davidson, 
1972), this plasma parameter was also examined. The smoothed 
average plasma ion density was available during the period 
1967 - July 1970 from Vela 4 , Pioneer VI and Vela 5 (Bame, 
Private communication). During June 19^7 - May I969 the data 
show a small positive correlation between n and the 60 MeV 
proton intensity; n decreases by 20^ while the proton intensity 
decreases by a factor of 2 . During the transition period 
there is a break in the correlation. However the absence of 
data during the recovery period does not allow us to comment 
on the reality of this effect. Further measurements of ion 
plasma density during the decreasing phase of the solar 
activity, magnetic field power spectra over the solar max- 
imum and the cosmic ray intensity gradients provided by 
Pioneer 10 , should soon be available and provide a more 
definite clue to the process of the cosmic ray modulation. 
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Flgure Captions 


Fig. 1 - Time history of the 6o MeV/nuc cosmic ray protons and 
alpha particles for the 1967-1971 period. The upper 
curve is the Deep River Neutron Monitor rate averaged 
over the same interval as the corresponding alpha and 
proton measurements . The vertical lines mark the 
maximum and minim-urn neutron monitor rates during a given 
interval. The 1965 60 MeV proton level is derived from 
Kinsey (1970). 

Fig. 2-60 MeV/nuc proton (a) and Alpha (b) data from Fig. 1 
plotted as a function of the corresponding Deep River 
Neutron Monitor rate. The dashed line is the theore- 
tical fit obtained for both model I and model II using 
the parameters listed in Table I. The indicated symbols 
divide the data into 5 different time periods . 
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Case II: Change in the effective size 
of the modulation region 

K(p,P,r,t) = pKo(P,r,t)P 

K^(P,r,t) = A(t) r < D(.P,t) 

Ko(P,r,t) = A(t)Exp(^) r > D(P,t) 

D(P,t) = 


A(t)xl0-21 
(cm^/ sec-GV) 

D(t) for 60 MeV 
Proton (A.U.) 

lo5-1.12 
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2 

0.9-1.12 

2-5 

1.12-1,8 

5-13 

1,8 - 2.7 
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